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ABSTRACT 

We have identified a candidate ~ 1 — 2 x 10 5 year old luminous white dwarf in NGC 1818, a young 
star cluster in the Large Magellanic Cloud. This discovery strongly constrains the boundary mass M c 
at which stars stop forming neutron stars and start forming white dwarfs, to M c <; 7.6 M . 
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1. INTRODUCTION 



Stara of mass comparable to that of the 3un evolve to 
form white dwarfs while, above some critical mass, Mr* 



stars detonate as type II supernovae instead, leaving neu- 
tron stars as remnants. The determination of M c is im- 
portant for understanding stellar evolution. It also has 
implications for the chemical evolution of galaxies, in that 
it affects the rate of enrichment of the interstellar medium 
and the total population of pulsars. The relationship be- 
tween the mass of a white dwarf and that of its progen- 
itor, as a function of the metallicity of the progenitor, is 
also important for understa n ding stellar evolution in the 
early universe (Jeffries 1997, van den Bcrgh fc Tammann| 
1991). Predictions of M c range from 6—10 M Q depend- 
ing o n the details of models and on t he metallicity of the 
star (Weidemann 199C. Jeffries 1997, Garcfa-Berro et al 



lOOT r Pols ct al., 1008 )t Almost no observational con ' 



strainto on M c are available except from one Galactic open 



cluster in which the m ost luminous of four whit e dwarfs 
implies M c > 5.7 M ( [Koester fc Reimers 1996| ). Statis- 
tical extrapolations of these da ta suggest M c ~ 8 M but 



with large formal uncertainty (Jeffries 1997). 

The young rich star clusters in the Large Magellanic 
Cloud (LMC) are particularly well suited for studies of 
the evolution of intermediate mass stars. They typically 
contain an order of magnitude more stars than Galactic 
open clusters, and membership determination is not gen- 
erally problematic as it is with Galactic clusters, which 
are often superposed against a dense curtain of disk stars. 
The clusters with age ~ 2-4x 10 7 yr have main -sequence 
turnoff masses ~ 7.5 — 9.5M , (Will et al., 1995) interest- 
ingly close to the theoretical range of values for M c . 

Models show that turnoff sta rs of this mass spend 
2xl0 6 yr on the red giant branch (Pols ct al., 1998 , SchaUea 
et al., 1992) before ejecting their hydrogen envelope and 



leaving rapidly fading white dwarfs (assuming they in fact 
have masses < M c ). Rich LMC clusters typically have 
~ 10 — 20 red giant branch stars, so we would expect 
the youngest remnant star to be ~ 1 — 2 x 10 5 years old. 
The models of Wood give luminosities for white dwarfs of 
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this mass an d age of 30 — 100L(T), w i th temperatures of 



100,000 K (Wood 1992, Wood 1991, D'Antona fc Mazz- 



itclli). Cooling is dominated by neutrino emission, and the 
remnants fade proportionally to ~ i~ 2 - 5 or faster, depend- 
ing on model details. 

Adopting a distance modulus for the LMC of (V — 
My) — 18.5 (Panagia et al., 1991) would imply apparent 
magnitudes for these young white dwarfs of V ~ 18 — 19.5. 
They would be visible even from the ground, although 
because of crowding, Hubble Space Telescope (HST) ob- 
servations would be required for accurate photometry. 
The white dwarfs would be distinguished from compa- 
rably bright main-sequence stars by their extreme blue- 
ness. To investigate the possibility of finding young white 
dwarfs in LMC clusters, we analysed HST archive images 
of NGC 1818. This cluster has mass ~ 2.8 x 10 4 M^, core 
radius r c = 2.0 pc, half-mass radius ~ 14 pc (Elson ct 
al., 1987| ), and age - 2 - 4 x 10 7 yr ( |Will et al., 1995| ). 
There are currently ~ 16 stars in the red giant phase. 

2. OBSERVATIONS 

The images of NGC 1818 were obtained with the Wide 
Field and Planetary Camera (WFPC2) on 1995 Decem- 
ber 8, with the F336W (= C/ 336 ), F555W (= V555), and 
F814W (= I S u) filters. Total e xposure times are 96 0, 880 
and 1290 seconds respectively ( Hunter et al., 1997| ). The 
images in each filter were coadded with a median filter to 
eliminate cosmic rays. 

DAOFIND was used to detect objects 4cr above the 
background. Point-spread function (PSF) fitting was used 
to help eliminate spurious detections: these include pri- 
marily structure in the PSF, particularly around saturated 
stars, and bright pixels along diffraction spikes. The fi- 
nal photometry was performed using an aperture with ra- 
dius 2 pixels, and aperture corrections and zero points 
were applied (Holtzmann ct al., 1995a| , Holtzmann et al 



1995b| , [Elson et al., 1998a| ). A value of E(B - V) = 0.05 
was adopted for the reddening. 

A color-magnitude diagram (CMD) in V555 vs {U^q — 
U555) for the Planetary Camera (PC) revealed a prominent 
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sequence of binary stars, and we explore the imp lications 
for th e bin ary population of NGC 1818 elsewhere ( Elson ct 
al., U 98a ). CMDs for the three WFC chips combined are 
shown in Figs, la and b. These cover a radial range from 
~ 2 core radii, out to about half way to the edge of the 
cluster. Stars with V555 < 17.5 are saturated; the main- 
sequence turnoff is at V555 « 14. (We note the presence 
of an apparent gap in the main-sequence at V555 ~ 20.3 
which may be due to a possible jump in stellar magni- 
tudes around ~ 1.5 — 2 due to the onset of convective 



overshooting ( Pols et al., 1998| ).) 

Models suggest that white d warfs should have (U — 
V) ~ -1.5 and (V-J)o ~ -0.4 ( [Wood 1991| , |Cheselka et] 
al., 1993); these values are roughly independent of age and 
metallicity over the range of interest. For white dwarfs the 
model colors are equiv alent to (V555 — Isia) — — 0.44 and 
(U 33 6 ~ V 555 ) = -1.85 flHoltzmann ct al., 1995b| Fig. 11). 

The boxes drawn in Figs, la and b delimit plausible 
ranges of colors and magnitudes for a young white dwarf. 
We inspected visually all detections within these boxes, 
and in a similar range of magnitudes and colors for the 
CMDs for the PC chip. All but one of the objects turned 
out to be spurious detections as described above. One, 
however, turned out to be unambiguously stellar. It is in- 
dicated with a filled circle in the CMDs, and is circled on 
the image of the cluster shown in Figure 2. 

Could the photometry, or identification as a white dwarf, 
be in error? Poisson error bars arc included in Fig. lb, and 
indicate that the extreme blueness of this object cannot be 
due to random photometric errors. Nor is it likely that the 
colors are biased bluewards by a residual cosmic ray in the 
F336W image: six independent images were coadded, so 
a cosmic ray would have had to hit the same pixel in sev- 
eral of them. Goodness-of-fit parameters (DAOPHOT's 
'sharpness' parameter and x 2 ) derived from fitting a PSF 
to this object are within the range expected for a star: 
there is no evidence that it is strongly peaked, as a cosmic 
ray would be, or resolved, as a background galaxy would 
be. (A background galaxy this bright would be easily iden- 
tifiable in a WFPC2 image.) 

The object is extremely unlikely to be a quasar: quasars 
of thi s magnitude have surface density ~ 1 per square de- 
gree (Collcss ct al., 1991, Boyle ct al., 1988j ) , so the proba- 
bility of finding one in our field is ~ 10~ J . It is even more 
unlikely to be a foreground star in the halo of our Galaxy: 
the number of Galactic stars of this magnitude expected 
in a WFP C2 field of view at high galactic latitude is less 
than one ( Santiago, Gilmorc fc Elson 1996] ). At a distance 
of, say, 10 kpc, a Galactic white dwarf would have abso- 
lute magnitude My ~ —3.5, which is far too bright to be 
an old white dwarf. On the other hand, given the num- 
ber of stars on the giant branch, we would expect to find 
about one young white dwarf if the turnoff mass in NGC 
1818 is indeed less than M c . Nor should we be surprised 
to find it outside the core of the cluster: more than half 
the red giants in the cluster are outside the core, and four 
are further from the cluster center than our white dwarf 
candidate. 

If we see one white dwarf with V ~ 18.5, would we ex- 
pect to see many more fainter ones? This depends on the 
cluster mass function and on the difference between the 
turnoff mass and M c . Given the number of red giants in 



the cluster, and assuming that M tU rnoff ~ M c , we might 
expect there to be 2 ± 2 other white dwarfs with magni- 
tude up to ~ 4 mag fainter than our candidate (ie. with 
18.5 < V < 22.5). Given that sample of stars in the PC 
image is incomplete due to the presence of many satu- 
rated stars, and that the WFC chips cover just over 50% 
of the outer parts of the cluster, the expected numbers 
are consistent with our observations. If M c is significantly 
greater than the turnoff mass of the cluster then stars with 
slightly smaller masses than the current turnoff mass will 
have evolved to form white dwarfs which by now will have 
cooled and faded. We might therefore expect there to be 
4-8 other white dwarfs at fainter magnitudes still; the ex- 
act number expected is dominated by small number statis- 
tics when the mass function is extrapolated. Due to the 
rapid fading, we would expect these to be ~ 5 magnitudes 
fainter (ie. V ~ 23.5) and they would be difficult or impos- 
sible to detect in the currently available WFPC2 images. 
The luminosity function of white dwar fs expected in young 
LMC clusters is discussed further in Elson et al., 1998a. 



Deeper images of NGC 1818 are scheduled during HST 
Cycle 7 (Project 7307) to look for a sequence of fainter, 
older white dwarfs. 



3. DISCUSSION 

What constraints on M c can we infer from the presence 
of a young white dwarf in NGC 1818? M c must be greater 
than the turnoff mass, so the main task is to determine an 
accurate value for this. The main uncertainties are not in 
determining the magnitude of the turnoff, but in convert- 
ing this magnitude to a mass. This requires a knowledge 
of both the metallicity and of whether stellar evolution 
models with or without convective overshooting are more 
appropriate. The metallicity of NGC 1818 has been de- 
tcrmincd from sp ectra of two stars to be [Fe/H]= —0.8 
(Will et al., 1995). This value is, however, lower than the 



expected value for the young population in the LMC (Ol- 
szewski et al., 1991), which is thought to be more like —0.2 
The membership of one of the stars is doubtful, and —0.8 
may in fact be too low. Two models with convective over- 
shooting, and metallicities [Fe/H]= —1.3 and —0.4 give 
an age for NGC 1818 of 4 x 10 7 yr, and a turnoff mass 
7.6±0. 1 M Q . Two models without overshooting giv e an age 
2.2 ±0.2 x 10 7 yr, and a turn off mass 9.0±0.5M o ( [Will ct 
al., 1994|Schaller ct al., 1992| ). From our own stellar evoul- 
tionary models, which include convective overshooting, we 
predict M c = 7.0 for [Fe/H]= -0.8 and M c = 7.7 M for 
[Fe/H]= —0.3. The convective overshoot parameter for 
these models is A c = 0.28 for M ~ 7.5 M©, which is close 
to but slightly higher than the models of Schaller et al.used 
by Will et al. Other groups use different values for A c , and 
a full discussion of the effect of the choice of thi s param 
eter on ou r results is deferred to a future paper (Elson c 
al., 1998a| ). Our observations imply M c <; 7.6 M©, ani 
perhaps M c <; 9.0 M©. A better determination of the 
metallicity of NGC 1818 would help constrain the value 
of M c further, as would deeper images in which any older 
white dwarfs would be visible. 

Another uncertainty affecting observational constraints 
on M c is the possiblity of an age spread among the stars 
near the main-sequence turnoff. For example, it is possi- 
ble that the lower mass stars formed first, by a sufficient 
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margin that some stars slightly below top of the main- 
sequence are evolving on to the red giant branch at the 
same time as more massive stars. The most massive stars 
have main-sequence lifetimes about 5 Myr, so the total age 
spread is at most a few Myr. Such a spread in age would 
mean that stars evolving on to the red giant branch could 
have masses ^ 0.3 M Q less than the measured turnoff mass 
which would imply M c <; 7.3 M Q . Measuring such an age 
spreads requires accurate photometry for stars just below 
the turnoff, and we will be acquiring this during Cycle 7. 
Another source of uncertainty is the possibility that the 
object is a member of a mass-transfer binary, but this is 
unlikely and would require it now be a tight white dwarf 
- neutron star binary. 

Because young white dwarfs are bright and very blue, 
detecting candidates even in ground based data is not dif- 
ficult. They would be expected to be present in any star 
forming region containing significant numbers of stars with 
masses near M c . For example, a recent CMD of an as- 
sociation in the Small Magellanic Cloud containing stars 
with ages 10 — 60 Myr contains three stars whose colors 
and ma gnitudes are consistent wi th those of young white 
dwarfs ( [Demers fc Battinelli 1998| ). 

Future HST observations of NGC 1818 and other young 
LMC clusters, to determine the cooling sequence, ages and 
masses of the white dwarf population in these clusters, 
should allow us to determine M c more precisely, and pos- 
sibly for a range of metallicities. 

4. SUMMARY 

We have identified a candidate luminous white dwarf in 
the young star cluster NGC 1818 in the LMC. The ob- 
ject is ~ 35 arcsec from the cluster center (about 4.5r c 
and 0.6r h ). It has coordinates 5:04:13.8, -66:26:33.4 
(J2000). In the HST passbands it has V 555 = 18.43, 
(V 55 5 - hu) = -0.26 and (U 336 - V 555 ) = -1-67. In the 



Johnson-Cousins system this corresponds to V = 18.44, 
(V - I) = -0.25 and (U - V) = -1.32. These values 
are corrected for reddening assuming E(B — V) — 0.05. 
Posisson errors are ±0.03 for V and ±0.04 for the colors. 
These do not include uncertainties in the transformation 
to the Johnson-Cousin system. The temperature is prob- 
ably ^ 20, 000 K but is poorly constrained by the (U — V) 
color. With the adopted distance modulus of 18.5, the 
object has absolute magnitude My = —0.06. 

If this object is indeed a white dwarf, th en its mass 
is probably 1.1 - 1.3 M Q ( |Weidcmann 1990[ ). The com- 
position of white dwarfs formed from high mass progeni- 
tors is expected to be Oxygen-Neon-Magnesium, but may 
be Carbon/Oxygen. A spectroscopic determination of its 
composition is a priority. If spectroscopic followup obser- 
vations confirm the identity of the candidate star as a lu- 
minous young white dwarf, we have strongly constrained 
the critical mass at which stars stop evolving to type II 
supernovae to M c <; 7.6 M Q . 

Note: A preliminary spectrum of the white dwarf candi- 
date was obtained at the AAT 5 March 1998. The spec- 
trum rules out the possibility that the object is a quasar. 
The velocity is indistinguishable from that of two other 
cluster members so it is also very unlikely to be a fore- 
ground object. Detailed modelling of the spectrum is cur- 
rently in progress. 
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Fig. 1. — Colour-magnitude diagrams for stars in NGC 1818 in the three WFC chips. Magnitudes are in the HST pass- 
bands and have been dereddened assuming E(B — V) = 0.05. Stars above V555 ~ 17.5 are saturated. The boxes indicate 
the range of colors and magnitudes where we might expect to find white dwarfs. All points in these boxes are spurious 
detections except the filled circle which is a candidate young white dwarf. Poisson errors in the {Uzzq — V555) colors are 
shown. Much of the scatter redwards from the main-sequence at V555 < 21 is due to the presence of binary stars. 



Fig. 2. — WFPC2 image of the inner regions of NGC 1818 in the F555W passband. The field is 130 arcseconds on a side. 
The core radius of the cluster is 8 arcsec. The candidate white dwarf is circled. (The diameter of the circle is 7 arcsec.) 
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